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such a metric. It turns out that the homogeneous models can only be either FLRW 
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well as hypersurface orthogonal spacelike Killing vectors) isometry. The latter can 
possess singularities of various kinds or none. The non-singular family is however 
unique and cylindrically symmetric. 
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1. Introduction 

It is to the credit of general relativity (GR) that the study of the Universe 
as a whole has become one of the most active areas of scientific research and it 
goes by the name relativistic cosmology. The models of the Universe we consider 
are given by the exact solutions of Einstein's field equations for gravitation.. The 
matter content of the Universe is taken as perfect fluid. The generally accepted 
model is the Friedman-Lemaitre- Robertson- Walker (FLRW) model which describes 
a homogeneous, isotropic and expanding Universe. Howsoever successful this model 
be, its homogeneity and isotropy are very special properties which could by no means 
be considered as generic enough for the Universe. It would hence be important to 
find cosmological solutions of Einstein's equations without imposing the conditions 
of homogeneity and isotropy. 

The first step in this direction came in the form of Bianchi models that 
are homogeneous but anisotropic. Inhomogeneity was the main concern of 
what is known as G2 cosmologies [1-3]. By G2 we mean a spacetime that 
admits a two parameter orthogonally transitive group of isometrics, that is, there 
exist two spacelike Killing vectors that are mutually as well as hypersurface 
orthogonal. In proper coordinates the metric depends upon only one space variable. 
Inhomogeneous models have been considered by some authors [4-7] of which the 
one due to Senovilla [6] is the most remarkable for being singularity free and yet 
having an acceptable physical behaviour. This was a startling result for it was 
generally believed on the strength of the powerful singularity theorems [8] that 
the occurrence of singularity is inescapable in GR so long as reasonable physical 
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conditions are satisfied. Here was apparently a counterexample to the theorems. 
Subsequent to the discovery it was found that the theorems became inapplicable 
because the solution in question did not obey one of the assumptions; the existence 
of causal compact trapped surfaces. Not only were all physical and and geometrical 
parameters are finite and regular for whole of the spacetime, the metric was shown 
to be geodesically complete [9] exhibiting the absence of a singularity of any kind. 

All prior attempts to manage the big-bang singularity or to construct a non- 
singular cosmological model involved either unphysical behaviour for matter, like 
p < 0, or quantum effects and new fields or modification of GR [10,11]. Senovilla's 
[6] was the first exact solution of Einstein's equations free of any kind of singularity 
and possessing all physically acceptable properties. Then the question arose, was 
it an isolated solution or did there exist a family of non-singular models? Ruiz 
and Senovilla [12] considered the general G2 metric separable in space and time 
in comoving coordinates and identified a large family of non-singular spacetimes 
with cylindrical symmetry. Is cylindrical symmetry necessary for a non-singular 
cosmological model? This is a pertinent question to be addressed next. 

In this paper we consider a general orthogonal metric separable in space and 
time in comoving coordinates and examine, in all generality, the permissible fluid 
models. We show that the already identified non-singular family [12] is unique. It 
turns out that the requirement of perfect fiuid imposes G^jGq and G2 isometrics 
on the spacetime; the first two are homogeneous Bianchi I and FLRW models while 
the last one alone can sustain inhomogeneity. It turns out that inhomogeneous 
models could be with or without singularity. However, the non-singular family is 
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unique and is cylindrically symmetric. We can thus characterise all the fluid models 
described by an orthogonal and separable metric [13]. 

In Sec 2 we set up the field equations for a perfect fluid for an orthogonal 
metric and in Sec 3 we prove a theorem characterising perfect fluid cosmological 
models and establishing uniqueness of the non-singular family. We conclude with a 
discussion. 

2. Field Equations 

We consider the orthogonal metric, 



ds^ = Ddt'^ - Adxj - Bdxl - Cdxl (2.1) 

with velocity field given by tt = VOdt. In comoving coordinates we assume the 
metric to be separable, ie ^ = A{xa)A{t) etc. The separability can invariantly 
be characterised by (i) = and (ii) a/d being constant over the 3- 

hypersurface. This can easily be verified from the following expressions for the 
kinematic parameters 



n_ 1 Mo Bo Co 

2Vd\A ^ B ^ C 



(2.2) 
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where = dA/dt and Do, = dD/dx"^. 

We note a general result arising out of the following two relations [14] , 



where 9,a,uj,Ua are the kinematic parameters; expansion, shear, rotation and 



We infer from the above relations: 

Lemma : In the absence of shear and vorticity, the expansion of fluid is constant 
over the 3-space orthogonal to the fluid congruence and, further, the acceleration 
also vanishes when the quantity g/goo is a separable function of space and time in 
comoving coordinates. 

Corollary : For the vorticity free spacetime with separability ( as is the case for the 
metric (2.1)), acceleration can be non-zero only if shear is non-zero. 

According to the Raychaudhuri equation [15], in the absence of vorticity 
acceleration is necessary for halting the collapse to avoid the singularity which in 
our case can only exist if shear is non-zero. Thus non-singular solutions represented 
by the metric (2.1) will always have to be both inhomogeneous and anisotropic. 

Einstein's field equations are 



(2.5) 




(2.6) 



acceleration, tta = Ua;iu'^- Note that we have assumed Ui = y/gooS^. 



Rik — -^Rgik = —^T^Tik, 
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(2.7) 



where for a perfect fluid 



Tik = (p + p)uiUk - pgik (2.8) 

The exphcit expressions for T^^ [16] look quite formidable and rather 
intimidating. Fortunately, we have discovered an underlying order in them that 
allows us to write the rest of them from a given two (one each of diagonal and off 
diagonal) by prescribing the appropriate permutation rules. We begin with 
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where a subscript denotes partial differentiation and here the assumption of 



separability is not effected. 



The successive cyclic permutations A — > S — > C — > ^ and 1 — > 2 — > 3 — > 1 will 
give T§,T§ from ; T|,Tf from and T^^T^ from T^^ To write from 
TqI, let 0^i2 (i.e. Aq iA2, ^ iT^) and B ^ C ^ D ^ B while follows 
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from for 2 ^ 3 ^ 1 ^ zO ^ -2{Tl ^ T^) and A ^ D ^ B ^ C ^ A. Thus 
we can write all ten T-^, given the two, one each of diagonal and off diagonal. 

3. Characterisation and Uniqueness 

The conditions implied by the perfect fluid character of the source are: T^q — 
0,Taf3 = ioT a ^ P and = Tl = . We shall a priori assume no isometries 
of any kind except the separability of the metric (2.1) in comoving coordinates. 
The implementation of the fluid conditions will lead to G3 (homogeneity) and Gq 
(both homogeneity and isotropy) symmetries for homogeneous models and G2 (only 
admitting two spacelike Killing vectors) symmetry for inhomogeneous models with 
or without singularity. Then all the fluid models described by the metric (2.1) are 
characterised and uniqueness of the already identified family of non-singular models 
is demonstrated. We prove the following theorem [13]. 

Theorem : The separable metric (2.1) can only represent the following kinds of 
perfect fluid cosmological models : 

(a) if homogeneous, then Bianchi I and FLRW models, 

(b) if inhomogeneous, then models with or without singularity. 
Further the non- singular family as already identified [12] is unique. 

Proof : Let us first of all write the three representative equations; T\q = 0, T12 = 



and Tl = T|, 
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/ No isometry : We assume no isometry to begin with. From eqn. (3.1), T^o = 0, 
wiU imply 



Di ^ D2 „ Ds ni 



(3.4) 



where 



Pi + kiQi = m, Bi = kiCi 



(3.5) 



and 



Pi = 



Co/C-An/A 



Qi 



Bo/B Ao/A 



Bo/B + Co/C Bo/B + Co/C 

others follow by the cyclic permertation. Here Ua and ka are constants. 

We integrate the exact differential 



(3.6) 



d{lnD) — {lnD)idxi + {lnD)2dx2 + {lnD)sdxs 



(3.7) 



along two different paths to give 



D{x^)^C^'{x^), A = C^/'=^(a;«) and 5 = C^^^ (a;«), 



(3.8) 



that is, the space dependence of the metric is all but determined. It remains to find 
C{xa)- Further for the time dependence we get 



and 77-1 = 1 + A;i + l/k2- 

It may be noted that we have so far used only the three equations Tao = 

to obtain the relations (3.8) and (3.9) which leave only C{xa) and two of 
A{t), B{t),C{t) to be determined. Substituting (3.8) in (3.2) and its permutants 
leads to C{xa) = const.. Thus the metric can only represent a homogeneous Bianchi 

1 model. 

On the other hand when Aq/A = Bq/B = Cq/C, the shear vanishes and so 
does the acceleration. The spacetime is then both homogeneous and isotropic which 
determine FLRW uniquely [17]. When Aq/A = Bq/B ^ Cq/C, eqn. (3.1) and its 
permutants will imply either Bianchi I or the spacetime admits a Gi isometry. This 
is the case we consider next. 

// Gi isometry : Let djdxj, be the spacelike Killing vector and hence the metric 
is a function of only two space variables, x\ and xi. 



+ /ci)^ + (1 + k^)^ + (1 + k^k^y^ = 



(3.9) 
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Note that eqn. (3.3) has the form 



where /i, /2, /s are functions of Xq,, containing respectively derivatives with respect 
to xi,X2 and 0:3. In this case /s = and T30 = 0. Eqn. (3.10) gives rise 
to two cases: (i) Aq/A = Bq/B ^ Cq/C and (ii) Aq/A ^ Bq/B ^ Cq/C and 
/i = const., /2 = const.. That means a is non-zero to give non-zero acceleration. 
It could be a viable case for a non-singular model as well. In (i) we can set 
^ = S for A{t) = B{t) is implied by Aq/A = Bq/B (a constant multiple can 
always be absorbed) and A[xa) = B{xa) can be done by an appropriate coordinate 
transformation. Eqn. (3.1) implies C{xa) — D^{xoi) and eqns. (3.1) - (3.3) give 
three equations to determine the space dependence of the metric. The Lie group 
analysis of the equations (see Appendix) leads to the inference that the functional 
dependence can only occur in the form A{x\ + X2).,A{x\ + X2) and A{xi/x2). 
The first two cases reduce to single variable dependence by suitable coordinate 
transformation, which we consider separately. The last case is obviously singular 
and could not be considered as a viable case for any kind of cosmology. 

In (ii) Aq/A ^ Bq/B ^ Cq/C, following the same route we get from 
2aO = 0; A, S and D in terms of C{xa) as before. Then T12 = determines C{xa) = 
ifi^i) + f{x2)Y- Eqn. (3.10) represents two equations; f\ = const., /2 = const. 
and two more similar equations. These will ultimately determine C{xa) = const. 
and again the spacetime is Bianchi I. 
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Thus Gi symmetry does not yield a viable fluid model. 

/// G2 isometry: Finally we have the spacetime general enough to sustain an 
inhomogeneous fluid. All fluid models are inhomogeneous and anisotropic. In view 
of eqns. (2.5) and (2.6), it follows that inhomogeneous spacetime has necessarily to 
be anisotropic. Inhomogeneous models can have singularities of different kinds or 
none. 

Ruiz and Senovilla [12] have thoroughly analysed this case and have shown that 
the spacetime possesses a rich singularity structure. The metric (2.1) will have 
models with singularity but not always of the big-bang kind as well as models free 
of singularity. The latter family is shown to be unique and cylindrically symmetric. 
Since non-singular solutions are allowed only in this case, the identifled non-singular 
family is unique for the general orthogonal metric (2.1). 

This completes the proof of the theorem. 

The most general non-singular metric [12] is given by. 



ds'^ = cosh^+^(at) cosh"-^ (nar) dt'^ - 



sinh^ (nar) 




p2 




cosh^ "(at) 
cosh " (nar) 



(3.11) 
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and K,n,a are constants. The coordinates range as — oo < t, z < oo, < r < 
oo, < < 27r and the metric has cylindrical symmetry. 

The fluid parameters are given by 



87rp = X 



(n - 1) (2n - 1) (n + 3)K cosh"^ (nar) + (n + 1) (n - 3) cosh~^ (at) 



(3.13) 



87rp = X 



(n- l)^(2n- l)i^cosh ^ (nar) + (n + l)(n - 3) cosh '^{at) 



, (3.14) 



where 



X = ^ cosh-(^+")(at) cosh'-"(nar). (3.15) 

Both density and pressure are positive and p < p for K >0. The equation of state 
p = 3p for radiation is admitted when n = 3 . Senovilla's model [6] further requires 
K = 1. The case K = gives the stiff fluid equation of state, p = p [18]. The case 
K = 1 has been considered separately and it has been shown that radial heat flux 
can be incorporated without disturbing the singularity-free character of the metric 
[18,19]. 

4. Discussion 

The main result of the paper is that the already identifled family of non- 
singular cosmological models is unique not only for the G2 metric but also for 
the general orthogonal metric separable in space and time in comoving coordinates. 
Thus the complete set of non-singular solutions has been identified. In the process of 
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establishing this result we have also been able to characterise all perfect fluid models 
described by the metric (2.1). They are: homogeneous Bianchi I and FLRW; and 
inhomogeneous with or without singularity (different kinds of singularities occur 
[12]). We assume no isometrics a priori, the perfect fluid conditions imply G3 and 
Gq symmetries for homogeneous, and G2 for inhomogeneous models. 

The non-singular character of fluid models singles out cylindrical symmetry. 
Like inhomogeneity, it is only a necessary condition but not sufficient. A kind of 
formal connection can be indicated between the non-singular metric (3.11) with 
K = 1 and the FLRW open model. The former can be thought of as arising out of 
a natural inhomogenisation of the latter [20] . The unfortunate feature of the metric 
(3.11) is that anisotropy does not decay with time (since a/d = const), which means 
it can never evolve into FLRW. There may, however, occur a non-singular solution 
when the assumption of separability is dropped, which may isotropise to FLRW at 
late times. That would be a very significant result for cosmology, but the situation 
becomes mathematically formidable. We are currently investigating this question 
for spherical and cylindrical symmetry. 
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Appendix 

An nth order system of differential equations 



E(x,y,y',...,yW) = 



(^.1) 



will have the Lie (point) symmetry [21] 



G = ^i{Xi,yi)-^ + r]i{Xi,yi) — 



{A.2) 



iff 



^["Ie = 0, 



(^.3) 



where G^"^^ is the nth extension of G needed to take care of the nth derivatives in 
(A.l) . 



We analyse the following three equations: 
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where a = (1 + 2A — A^)/(l + A) and I and A are constants. 
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The standard Lie analysis [22] gives the Lie point symmetries of the above 
equations as 



G2 = 4- (^-8) 

0x2 



To reduce the system of partial differential equations (A.4)-(A.6) to one of ordinary 
differential equations (and hence solve them) we need to decide on an appropriate 
independent variable. From (A. 7) and (A. 8) we have the possiblity of independent 
variables 

u = ClX\ + C2X2 (^-11) 

and from (A. 9) 

«=^. (A12) 

Xi 

The fourth symmetry implies that u = u{xi,X2)- 
Now, (A. 11) implies 

A = A{xi+X2) D^D{xi+X2) (A.13) 
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while (A. 12) implies 



A=-a(-) (AM) 

Xi \XiJ 

with 

B = .,p(|) (A15) 
coming from the addition of (A. 10). On the other hand, u = u{xi,X2) implies 



A = A{u), D = D{u) {A.I&) 

which further implies 



u = f{xi+xl). (A17) 

We could take other combinations of the symmetries. The only one of relevance 
is (A.9) + k (A.IO). This gives 

and is equivalent to (A. 14) and (A. 15) for the purposes of determining whether A 
and D contain singularities. 
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